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ABSTRACT 

The oxidation of DNA resulting from reactive oxygen 
species generated during aerobic respiration is a 
major cause of genetic damage that, if not 
repaired, can lead to mutations and potentially an 
increase in the incidence of cancer and aging. A 
major oxidation product generated in cells is 
8-oxoguanine (oxoG), which is removed from the nu- 
cleotide pool by the enzymatic hydrolysis of 
8-oxo-2 -deoxyguanosine triphosphate and from 
genomic DNA by 8-oxoguanine-DNA glycosylase. 
Finding and repairing oxoG in the midst of a large 
excess of unmodified DNA requires a combination 
of rapid scanning of the DNA for the lesion followed 
by specific excision of the damaged base. The repair 
of oxoG involves flipping the lesion out of the DNA 
stack and into the active site of the 8-oxoguanine- 
DNA glycosylase. This would suggest that thermo- 
dynamic stability, in terms of the rate for local de- 
naturation, could play a role in lesion recognition. 
While prior X-ray crystal and NMR structures 
show that DNA with oxoG lesions appears virtually 
identical to the corresponding unmodified duplex, 
thermodynamic studies indicate that oxoG has a 
destabilizing influence. Our studies show that 
oxoG destabilizes DNA (AAG of 2-8kcalmol 1 over 
a 16-1 16 mM NaCI range) due to a significant reduc- 
tion in the enthalpy term. The presence of oxoG has 
a profound effect on the level and nature of DNA 
hydration indicating that the environment around 
an oxoG«C is fundamentally different than that 
found at G«C. The temperature-dependent imino 
proton NMR spectrum of oxoG modified DNA 



confirms the destabilization of the oxoG»C pairing 
and those base pairs that are 5' of the lesion. The 
instability of the oxoG modification is attributed to 
changes in the hydrophilicity of the base and its 
impact on major groove cation binding. 

INTRODUCTION 

The oxidation of guanine to 8-oxoguanine (oxoG) in 
DNA is a frequent transformation that occurs in cells as 
a result of the metabolic generation of reactive oxygen 
species or exposure to agents that induce oxidative stress 
(1-5). It has been demonstrated that this modification is 
both cytotoxic and mutagenic (6-9), the latter presumably 
due to misincorporation errors by either replicative poly- 
merases (10) or by translesion synthesis (TLS) polymer- 
ases that can bypass lesions, albeit with reduced fidelity 
(11-15). As is frequently the case for the more common 
types of DNA lesions, there are specific proteins involved 
in the repair of oxoG, as well as those that restrict its 
incorporation into DNA via the dNTP pool. The latter 
process is performed by MutT that hydrolyzes 
2'-deoxy-oxoGTP to 2'-deoxy-oxoGMP and prevents its 
use as a substrate by DNA polymerases (16-18). The 
major repair pathway for removal of oxoG from DNA 
is base excision repair (BER) and the initial step is per- 
formed by 8-oxoguanine-DNA glycosylase (MutM in 
bacteria or OGG in mammals), which catalyzes the hy- 
drolysis of the lesion from the DNA backbone to afford 
the free oxoG base and leaving an abasic site in the DNA 
(19-21). To accomplish this reaction, the oxoG nucleotide 
must flip out of the DNA stack into the active site of the 
protein (22). After excision of the oxoG, the remaining 
abasic site is then serially processed by a number of 
other BER proteins to eventually restore a normal G«C 
base pair (23,24). 
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The mechanism by which the OGG glycosylase finds its 
substrate in the midst of a huge excess of undamaged 
DNA is an area of interest that has potential clinical im- 
plications since oxidative damage is associated with both 
aging (25,26) and cancer (27-29). Structural analyses of 
DNA with an oxoG»C base pair have been reported using 
a number of experimental methods. The crystal structure 
shows that the base pairing is very similar to a canonical 
Watson-Crick G»C base pair (30). Low temperature (5°C) 
NMR studies of a 5'-d(CCACTA-oxoG-TCACC) 
sequence also indicate a normal Watson-Crick B-form 
arrangement (31). However, the NMR studies provided 
evidence for a 3-5° C decrease in thermal stability (T M ) 
as measured by the temperature-dependence of the 
chemical shifts for the non-exchangeable protons, as well 
as for line-widths of the exchangeable imino protons. The 
origin of this destabilization was attributed to a local 
change around the oxoG»C base pair without a large-scale 
change in the overall structure. UV melting and DSC 
studies of 5'-d(GCGTAC-oxoG-CATGCG)-3'«3'-d(CG 
CATGCGTACGQ-5' duplex yielded the same general 
observation of a reduced (2-3°C) thermal stability in 
1 M salt concentration. Note that the NMR studies were 
done at 200 mM salt concentration and at higher strand 
concentration (32). The base pairing of 8-oxo-dG with dC, 
dT, dG and dA has also been measured by NMR, and it 
was concluded, consistent with both the high-resolution 
crystal and NMR structures, that the 8-oxo-dG»dC has 
a normal high-affinity Watson-Crick base pairing ar- 
rangement (33). 

Based upon structural studies that indicate normal 
Watson-Crick oxoG»C base pairing, it was proposed 
that OGG (and MutM) can detect its cognate lesion 
through 'recognition' of the oxoG lesion while it is 
stacked within the double helix (22,34-37). In this 
regard, the discrimination between oxoG«C and G»C 
could involve the unique H-bond donor available at 
N7-H on oxoG with an H-bond acceptor on the 
glycosylase (e.g. the amide carbonyl of Gly42). In 
summary, a series of discrimination events induce the 
specific helical extrusion and entry of oxoG into the 
enzyme's activie site (34). The proposed direct recognition 
of the intrahelical oxoG by OGG and MutM contrasts 
with the suggested mechanism by which uracil 
DNA-glycosylase (UNG) finds and repairs deoxyuridine 
(dU) in DNA (38-41). In the case of dU, the experimental 
data are consistent with UNG passively finding and 
trapping extrahelical dU residues rather than facilitating 
the selective extrusion of the cognate dU from the DNA 
stack base. Once inside the active site of their respective 
glycosylases, both oxoG and dU are efficiently hydrolyzed 
from the DNA leaving an abasic site. It is of course quite 
possible, that the two glycosylases use different mechan- 
isms to find their cognate lesions. 

The normal appearance of an oxoG»C base pair by 
crystallography and NMR is remarkable since there are 
major differences in the atoms that line the floor of the 
major groove at the oxidized lesion. Specifically, the elec- 
tronegative 7 -N 8 ~ atom and unpolarized C8-H on G are 
replaced on oxoG by a electropositive N-H 5+ and an elec- 
tronegative sp 2 -O b ~, respectively. Accordingly, we 



anticipated that the arrangement of cations and waters 
of hydration in the major groove near an oxoG»C base 
pair would be quantitatively and qualitatively different 
than at a canonical G»C, and that this would alter local 
DNA stability and its conformational dynamics. 
Previously, we observed that 7-deazaguanine caused an 
enthalpy-derived destabilization of DNA, which we 
attributed to the elimination of a major groove cation 
binding site (42,43). 

Because high-resolution structural studies generally do 
not account for water and cations, we initiated an inves- 
tigation of the thermodynamic properties of an oxoG»C 
base pair with specific interest in the effect on cation 
binding and hydration. To complement the thermo- 
dynamic studies, the temperature-dependent changes in 
the imino proton NMR spectrum were determined. The 
results show that there are significant thermodynamic dif- 
ferences between the unmodified and oxidized base, and 
that there is a reduction in base pair stacking and 
enhanced exchange of the imino protons at the oxoG»C 
base pair. The origins of this destabilization are discussed. 

MATERIALS AND METHODS 

Materials 

The unmodified and modified oligodeoxynucleotides 
were synthesized by the Midland Certified Reagent Co. 
(Midland, TX, USA) and purified by anion-exchange 
HPLC. Purity of the oligodeoxynucleotides (Figure 1) 
were verified by HPLC using a semi-preparative reverse- 
phase column (Phenomenex, phenyl-hexyl, 5 um, 250 x 
10.0 mm) equilibrated with 10 mM triethylammonium 
acetate (pH 7.0). All the oligomers were desalted on a 
G-25 Sephadex column, lyophilized to dryness and 
characterized by MALDI-TOF-MS (calculated mass for 
[M-H]- m/z 3662.4; found m/z 3663.88). The dry oligo- 
mers were annealed by dissolving the single-stranded 
oligodeoxynucleotides in appropriate buffer, heating the 
solution up to 85°C for 15min and allowing it to slowly 
cool down to room temperature. The oligomer concentra- 
tions for the NMR studies were determined by UV using 
an extinction coefficient of 1.11 x 10~ 5 M _1 cm _1 at 
260 nm and 25° C. The strand concentrations of the oligo- 
nucleotides (ODNs) for the thermodynamic studies were 
determined at 260 nm using the following molar extinction 
coefficients assuming similar extinction coefficients for 
oxoG and G (in mM-'cm -1 ): 95 (GCGAATTCGC), 



ODN Sequence 

1 5'-d(GCGAATTCGC)-3' 

2 5'-d(GCGAATTC-oxoG-C)-3' 

3 5"-d(GAGAGCGCTCTC)-3' 

4 5'-d(GAGA-oxoG-CGCTCTC)-3' 

5 5'-d(CGCGTTTTCGCG)-3' 

6 5'-d(CGCGTTTTC-oxoG-CG)-3" 

7 5'-d(GCGAATTC-c 7 G-C)-3' 

8 5"-d(GAGA-c 7 G-CGCTCTC)-3' 

Figure 1. Sequence of the ODNs; oxoG, 8-oxoguanine; c 7 G, 
7-deazaguanine. 
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105 (CGCGTTTTCGCG) and 112 (GAGAGCGCT 
CTC). These values were calculated by extrapolation of 
the tabulated values of the dimers and monomer bases 
(44) at 25 and 80°C using procedures reported earlier (45). 

Temperature-dependent melts by UV spectroscopy 

The UV melting experiments were performed on a Cary 
UV 300 Bio UV-visible spectrophotometer (Varian Inc., 
Palo Alto, CA, USA) equipped with a Cary Temperature 
Controller. The melting of the ODNs were monitored at 
260 and/or 275 nm in the temperature range of 0-90°C at a 
heating rate of 1.0°C/min with a data interval of 0.4°C. 
The transition temperature (T M ) was determined by the 
analysis of the first derivative and shape (46) of the 
melting curves. The nature of complex formation of the 
ODNs was studied by performing the UV melts as a 
function of their total strand concentration (3-330 uM). 
Differential binding of the counterions (A« Na+ ) and 
water (A« w ) molecules were determined by performing 
the UV melts of the ODNs in the presence of salt 
(10-200 mM NaCl), and ethylene glycol (0.5-3.0 m at 
lOmM NaCl), respectively, to understand the molecular 
changes accompanying the helix to coil transitions. These 
parameters are calculated using the procedure described 
earlier (47^19). 

Thermal denaturation by differential scanning calorimetry 

Calorimetric enthalpies for the helix to coil transitions of 
the ODNs were acquired on a VP-DSC differential 
scanning calorimeter (MicroCal, LLC, Northampton, 
MA, USA). In a typical DSC experiment, ~200uM of 
total strand concentration of the ODN solution in the 
sample cell (0.506 ml cell volume) was scanned from 
0-90°C at a rate of 45°C/h with buffer in the reference 
cell. A buffer versus buffer scan was also done under 
similar conditions and subtracted from the subsequent 
experimental runs. A 16 s filter period was used for all 
the data acquisitions. The resulting thermograms were 
plotted as heat capacity (AC P ) versus temperature 
profiles using Origin 7.0 software provided with the instru- 
ment. Analysis of the integrated plots of the anomalous 
AC P versus temperature curves (f AC p dT) and normaliza- 
tion for the number of moles, yields the molar calorimetric 
enthalpy (A// cai ). The molar calorimetric entropy (AS cal ) 
was calculated from integration of the AC P /T versus T 
curves [J(AC p /T)dT\. We use the Gibbs equation 
(Equation 1), or Equation (2), to calculate the Gibbs 
free energy at 20° C: 



AGcai(7) = Aflcai — TAScal, 



AGcai(7) = AHca 



T 



(1) 



(2) 



Pressure perturbation calorimetry 

A VP-DSC differential scanning calorimeter equipped 
with a pressure perturbation accessory was used to 
measure the heat (A0 resulting from the application of 
an additional pressure change (AP) above the solutions of 



the calorimetric cells. A complete description of this 
pressure perturbation calorimetry (PPC) technique can 
be found elsewhere (50-52). Prior to the PPC scans, 
standard DSC curves were obtained for the oligonucleo- 
tides using a 10 mM sodium phosphate buffer at pH 7, 
with and without the oxoG lesion, to determine the tem- 
perature range and temperature step to be used in the PPC 
experiment. A sample solution with 0.7 mM (control) or 
1 . 1 mM (modified) concentration in strands is allowed to 
equilibrate against the same buffer solution (reference cell) 
at constant temperature and external pressure. The 
external pressure is then increased by ~50psi, causing 
heat to be absorbed differentially by the sample and ref- 
erence cells. These heats, AQ, are obtained from integra- 
tion of the compression and decompression peaks 
resulting from switching the external pressure on and off 
at particular temperatures determined from the DSC 
curves. In addition, buffer-buffer, buffer-water and 
water-water scan are performed for the proper setting of 
the pre- and post-transitional baselines, and baseline sub- 
traction. The AQ values are used to measure the apparent 
coefficient of thermal expansion, a(7"), from integration of 
the relationship: 



OQ/dPh 



-T(dV/dT) P = -Ta(T)V 



yielding: AQ = — TV<x(T)AP; where V is the apparent 
molar volume of the solute. 
For a two component system: 

a(T) = a 0 -[AQ/(TVAP)], 

where a Q is the thermal coefficient of the solvent. 
Integration of a(7) over the temperature range of the un- 
folding reaction, Ja{T)dT, yields the relative volume 
changes of the solute, A V/ V, where A V is the unfolding 
volume of the macromolecule. The value of V is obtained 
from the following Equation (53): 

V= M/p 0 -(p 0 - p)/p 0 C, 

where M is the molecular mass of DNA, p a and p are the 
densities of the solvent and solution, respectively. The 
density of these solutions are measured with an Anton 
Paar (Graz, Austria) DMA densitometer in the differen- 
tial mode, using two 602-M micro cells, each with a 
volume of ~ 1 50 ul. The reference cell is filled with water 
while the measuring cell is filled with solution or buffer. 
The density, p, is calculated from the oscillation period T 
of the cell using the following relationship: 

p = AT 2 +B, 

where A and B are constants determined from calibrating 
densities (and periods) of water and air (53,54). 



NMR 

Samples for the observation of exchangeable protons were 
dissolved to a duplex concentration of 0.25 mM in 180 ul 
of 10 mM NaH 2 P0 4 , 200 mM NaCl, 50 uM Na 2 EDTA 
buffer (pH 7.0) containing 9:1 H 2 0/D 2 0 (v/v). One- and 
two-dimensional (ID and 2D) NMR experiments were 
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performed on a Bruker Avance spectrometer operating at 
600 MHz. Chemical shifts were referenced to water reson- 
ance. NMR data were processed using TOPSPIN software 
(2.0.b.6, Bruker Inc., Karlsruhe, Germany). One- 
dimensional NMR spectra for the exchangeable protons 
were recorded at 5, 15, 25, 35, 45, 55 and 65°C. The 'H-'H 
NOESY spectra of unmodified and modified samples in 
H 2 0 were collected at 5°C, with 70 and 250 ms mixing 
time and relaxation delay of 2.0 s (55,56). These experi- 
ments were recorded using a field gradient Watergate 
pulse sequence for water suppression (57). 

RESULTS 

Effect of oxoG on the transition temperature of the 
oligonucleotides monitored by UV spectroscopy 

Table 1 summarizes the thermodynamic parameters 
obtained from the UV melting analysis of the unmodified 
ODNs and the corresponding ODNs with an oxoG modi- 
fication. ODN-1 melts as a single transition from duplex 
to random coil with a T M of 29.5 and 53.0°C at 10 
(Figure 2A) and lOOmM salt (Supplementary Data) con- 
centrations, respectively. Replacing the G-9 with an oxoG 
caused a destabilizing effect and the Ar M relative to 
control ODN-2 was 6.3 and 13.8°C at 10 (Figure 2A) 
and lOOmM (data not shown) salt concentrations, re- 
spectively. In concordance with a duplex structure, the 
7m values of ODN-1 and ODN-2 showed strand depend- 
ence in the range of ~8-256 uM total strand concentration 
in the presence of 10 (Figure 3) and lOOmM 
(Supplementary Data) NaCl concentration. 

The duplex sequence ODN-3 (5'-GAGAGCGCTCTC) 
showed a T M value of 41.0 and 59.5°C in 10 (Figure 2B) 
and 100 mM NaCl, respectively. Modification of the G-5 
residue (ODN-4) resulted in a Tf^ reduction of 1.6 and 
8.3°C in 10 (Figure 2B) and lOOmM salt, respectively. 
Strand-dependent experiments done at lOmM salt con- 
centrations showed an increasing T M with an increase in 
the total strand concentration in the range of ~4-330 uM 
(Figure 3). A similar trend is observed at lOOmM salt 
concentration (Supplementary Data). This strand 



dependent single transition is evidence of a duplex to 
random coil denaturation. 

The 5'-CGCGTTTTCGCG sequence (ODN-5) forms a 
hairpin with a T 4 loop and shows a monophasic transition 
for random coil formation. The T M is 68.4 and 69.9°C in 
the presence of 10 (Figure 2C) and 100 mM 
(Supplementary Data) salt, respectively. The oxoG 
modified ODN-6 at the G-10 melted at 56.5 and 61.0°C 
in 10 (Figure 2C) and lOOmM (data not shown) salt, re- 
spectively. The strand dependence of the sequences 
(~3-264uM total strand concentration) done at 10 
(Figure 3) and lOOmM salt (data not shown) concentra- 
tions did not show any change in the T M as expected for 
the unfolding of an intramolecular hairpin. 

Comprehensively, it can be noted that, oxoG causes a 
T M reduction in all the studied ODNs, irrespective of the 
flanking sequences and the salt concentration. In addition 
to the change in T M , the hyperchromicity, which is an 
indicator of reduced base stacking, for the melts of the 
oxoG substituted ODNs is significantly reduced relative 
to the natural sequences. The drop in hyperchromicity 
for both ODN-2 and ODN-4 relative to ODN-1 and 
ODN-3, respectively, is ~40%, while the change is 
~30% for the hairpin structure. 

Calorimetric analysis of the oxoG modified 
oligonucleotides 

Representative DSC thermograms of the unmodified and 
oxoG modified ODNs in the presence of lOmM NaCl 
concentration are shown in Figure 4 and the correspond- 
ing thermodynamic parameters at 10 and lOOmM NaCl 
are given in Table 1 . The reported thermodynamic values 
represent an average of the analysis of five or six consecu- 
tive scans. 

Sequences ODN-1, ODN-2, ODN-3 and ODN-4 unfold 
as a single monophasic transitions as expected for a 
duplex to random coil formation. Both modified oligo- 
mers, ODN-2 and ODN-4, show a decrease in the T M 
compared to the corresponding unmodified ODNs at 10 
(Figure 4A and B) and 100 mM salt concentrations 
(Table 1) (Supplementary Data). The calorimetric 
enthalpy of ODN-2 yielded an endothermic difference of 



Table 1. Standard thermodynamic profiles for the formation of DNA duplexes at 20°C a 



ODN 


NaCl 






AG" 20 


TAS C , X 


A» Na + 


A«„ 




(mM) 


(°C) 


(kcal/mol) 


(kcal/mol) 


(kcal/mol) 


(moP 1 DNA) 


(/mol DNA) 


1 


10 


29.5 


-80.1 


-5.6 


-74.5 


-2.2 ± 0.2 


-30 ± 4 




100 


53.0 


-72.3 


-8.2 


-64.1 


-1.7 ± 0.1 


-22 ± 3 


2 


10 


23.2 


-44.4 


-2.3 


-42.1 


-1.7 ± 0.1 


-15 ± 1 




100 


39.2 


-33.2 


-3.0 


-30.1 


-1.2 ± 0.1 


-10 ± 1 


3 


10 


41.0 


-78.2 


-6.9 


-71.3 


-3.4 ± 0.2 


-41 ± 3 




100 


59.5 


-92.0 


-12.5 


-79.5 


-3.6 ± 0.2 


-43 ± 4 


4 


10 


39.4 


-39.3 


-3.2 


-36.1 


-1.0 ± 0.1 


-14 ± 2 




100 


51.2 


-39.5 


-4.6 


-35.0 


-0.9 ± 0.1 


-13 ± 2 


5 


10 


68.4 


-31.0 


-4.4 


-26.6 


-0.3 ± 0.1 


-18 ± 2 




100 


69.9 


-27.7 


-4.1 


-23.6 


-0.2 ± 0.1 


-16 ± 1 


6 


10 


56.5 


-16.6 


-1.8 


-14.8 


-0.3 ± 0.1 


-7 ± 1 




100 


61.0 


-7.8 


-1.0 


-6.9 


-0.1 ± 0.1 


-3 ± 1 



'All parameters are measured from UV (T M ) and DSC melting curves in 10 mM sodium phosphate buffer (pH 7.0) using ~10uM DNA. The 
observed standard deviations are: T M (±0.7), AH c . it (±3), AG° 2 o (±0.15), TAS ca \ (±2). 
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100 



Temperature / "C 



Figure 2. Absorbance versus temperature curves in 10 mM sodium 
phosphate buffer (pH 7) for — lOuM total strand concentration of 
the ODNs: (A) ODN-1 (open square) and ODN-2 (filled square); 
(B) ODN-3 (open triangle) and ODN-4 (filled triangle) at 260 nm; 
and (C) ODN-5 (open circle) and ODN-6 (filled circle) at 275 nm. 
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Figure 4. Calorimetric transition curves for ODNs in 10 mM sodium 
phosphate buffer (pH 7): (A) ODN-1 (open square) at 200 uM and 
ODN-2 (filled square) at 256 uM. (B) ODN-3 (open circle) at HOuM 
and ODN-4 (filled triangle) at 130uM. (C) ODN-5 (open circle) at 
200 uM and ODN-6 (filled circle) at 264 uM total strand concentration. 



CO 

o 




Figure 3. T M dependence on the strand concentration of ODNs in 
10 mM sodium phosphate buffer (pH 7): ODN-1 (open square) 
-8-178 uM, ODN-2 (filled square) -10-256 uM, ODN-3 (open 
triangle) -4-330 uM, ODN-4 (filled triangle) -4-127 uM at 260 nm, 
ODN-5 (open circle) -7-175 uM and ODN-6 (filled circle) 
—3-264 uM total strand concentration at 275 nm. 



35.7 and 39.1kcal/mol compared to ODN-1 at 10 and 
100 mM NaCl, respectively. Similarly the AA// cal values 
for ODN-4 was calculated to be 38.9 and 52.5 kcal/mol at 
10 and lOOmM salt concentration, respectively, compared 
to the unmodified ODN-3 at low and high salt. 

The oxoG modified hairpin sequence, ODN-6 also melts 
as a single transition at 10 (Figure 4C) and lOOmM 
(Supplementary Data) NaCl concentrations with the T M 
values as 56.5 and 61.0°C, respectively. The values are ap- 
proximately the same as obtained from the UV curves con- 
firming the structure of the ODN as a hairpin. The AA// ca i 
values obtained for the oxoG modified hairpin sequence at 
10 and lOOmM NaCl are 14.4 and 19.9kcal/mol, respect- 
ively. The lower AAH C£L \ values correspond to a single 
modification in this hairpin molecule while the duplexes 
have two modifications per complex. 

A decrease in the thermodynamic parameters for all the 
oxoG modified ODNs is reflected in the thermodynamic 
stability of the transitions represented by the Gibbs free 
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-1.5 -1.0 
log [Na + ] 

Figure 5. T M dependence on the salt concentration of ODNs (~7 uM 
total strand concentration) in 10 mM sodium phosphate buffer (pH 7) 
at 260 nm for ODN-1 (open square), ODN-2 (filled square), ODN-3 
(open triangle), and ODN-4 (filled triangle), and at 275 nm for ODN-5 
(open circle) and ODN-6 (filled circle). 



energy changes (AAG° at 20°C) of ~3.2 ± 0.6kcal/mol in 
10 mM NaCl (Table 1). This reduction in all the thermo- 
dynamic parameters for the oxoG modified ODNs 
compared to the corresponding unmodified ODNs 
confirms the UV melting results and the previous NMR 
and thermodynamic reports (31,32) that the introduction 
of oxoG in these sequences destabilizes DNA. The differ- 
ences in the magnitudes of the ODNs correspond to se- 
quences differences, i.e. the oxoG«C base pair is placed in 
different base-pair environments, 5'-C-oxoG-C and 5'-A- 
oxoG-C. 

Effect of counterions and water in the stability of oxoG 
modified oligonucleotides 

In most cases, the melting of DNA is a phenomenon 
accompanied with simultaneous release of water molecules 
and cationic counterions, and the latter can be hydrated. 
These counterions and water molecules are known to add 
to the stability of native DNA. Therefore, studying the 
hydration and cations associated with the different 
duplex structures provides information about the factors 
that contribute to the enthalpy and entropy terms respon- 
sible to the overall stability of the duplex versus random 
coil structures and for modified versus unmodified DNA. 
The melting of the oxoG modified ODNs (ODN-2, 
ODN-4 and ODN-6) was studied in the presence of 
10-200 mM NaCl to understand the effect of oxoG on 
the release/uptake of the Na + in these ODNs. 

ODN-2, ODN-4 and ODN-6 show a 17.0, 19.6 and 
9.0°C increase in thermal stability, respectively, upon 
raising the salt concentration from 10 to 200 mM. The 
linear dependence of T M on Na + concentration is graph- 
ically represented in Figure 5. Each plot gives a linear 
regression coefficient of ~0.99. The slope of these lines 
along with the A// cal and T M determined by calorimetry 
at 10 and 100 mM NaCl was used to calculate the change 
in the number of moles of Na + ion bound (A« Na+ ) upon 
the melting of the three ODNs (Table 1). At a salt con- 
centration of 10 mM, the values of An Na+ for the forma- 
tion of ODN-2, ODN-4 and ODN-6 are found to be -1.7, 



— 1.0 and —0.3, respectively (Table 1). Compared to 
ODN-1 and ODN-3, the values are lower with a difference 
of 0.5 and 2.4 in the case of ODN-2 and ODN-4, respect- 
ively, in lOmM salt, indicating lower counterion binding 
in the presence of oxoG in the DNA. The oxoG modified 
hairpin, ODN-6, did not show any change in the number 
of cations released within the experimental error (Figure 5 
and Table 1). Increasing the bulk salt concentration to 
100 mM led to a decrease in the A« Na+ values, which 
could be attributed to the higher screening of the phos- 
phates at a higher salt concentration (Table 1). 

The changes in the hydrophobic/structural water 
associated with the melting process were estimated by 
the osmotic stress method using ethylene glycol in the 
range of 0.5-3.0 m at 10 mM salt concentration. Spink 
and Chaires (58) showed that low-molecular weight 
co-solutes, such as ethylene glycol, influence DNA 
melting only through changes in water activity. Hence, 
the slope of a plot of T M versus logarithm of the activity 
of water (log a w ) allows the calculation of A« w , the change 
in the number of moles of water per base pair upon 
melting of the duplex. The thermal behavior of the 
ODNs in the presence of osmolyte in 10 mM NaCl is pre- 
sented in Figure 6. It is observed that with an increase in 
osmolyte concentration from 0.5 to 3.0 m, the T M shows a 
decreasing trend. A plot of T M versus log a w (Figure 7) 
gives negative slopes, each with a linear regression coeffi- 
cient of ~0.99. The number of water molecules associated 
with the transition of these ODNs was calculated using the 
values of A// ca i and T M obtained from the DSC curves at 
10 mM NaCl for each of the negative slope values 
obtained in the presence of ethylene glycol at this salt 
concentration. The values calculated are given in 
Table 1. The AA« W values obtained for the UV melt of 
ODN-2, ODN-4 and ODN-6 in the presence of ethylene 
glycol in lOmM NaCl are —15, —27 and —11, respectively. 
At the higher salt 100 mM concentration, the A« w values 
are lower (Table 1) due to higher screening of the water 
dipoles. 

In summary, the unfolding of each ODN is 
accompanied by a net release of counterions and water 
molecules. However, the extent of release is much lower 
with the oxoG modified ODNs. 

Effect of OxoG in the unfolding volume of the ODN 

PPC was used to determine the unfolding volume of 
ODN-4 and its unmodified control ODN-3 in order to 
obtain additional insight on the hydration effects accom- 
panying the substitution of an oxoG»C base pair for G«C. 
The PPC curves of ODN-3 and ODN-4 are shown in 
Figure 8. In this experiment, the temperature dependence 
of the coefficient of thermal expansion (a) is followed as a 
function of temperature. The magnitudes of a and of the 
slopes of the pre- and post-transition baselines (Table 2) 
indicate the type of chemical groups that are exposed to 
the solvent throughout the transition. At low tempera- 
tures, the positive values of a and the negative slopes of 
the pre-transition baselines, indicate the exposure of 
hydrophilic groups (i.e. sugar-phosphate backbone) to 
the solvent in their folded states. However, ODN-3 has 



Nucleic Acids Research, 2011, Vol. 39, No. 15 6795 



a lower a-value coupled with a less negative slope relative 
to ODN-4. This suggests that the duplex state of ODN-4 is 
slightly more hydrophilic than ODN-3, despite the fact 
that both duplexes have similar molar volumes (<t>^) 
(Table 2). At higher temperatures, due to the additional 
exposure of aromatic groups to the solvent, the 
post-transition baselines have lower a-values and slopes 



closer to zero; their magnitudes are more or less similar. 
This indicates, as expected, similar random coil states for 
both ODN-3 and ODN-4. 

The helix-coil transitions for each ODN are shown as 
negative peaks with T M s (Figure 8) similar to the T M s of 
the DSC unfolding curves. Integration and analysis (52) 
of these peaks yielded unfolding volumes (AV V ) of 
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Figure 6. UV melting curves for ODN-2 (A), ODN-4 (B) and ODN-6 
(C) in the presence of 0.5 (filled square), 0.85 (open square), 1.4 (filled 
circle), 2.4 (open circle) and 3.0 (filled triangle) molal ethylene glycol in 
10 mM sodium phosphate buffer, pH 7.0. 
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Figure 7. T M dependence on the ethylene glycol concentration of 
ODNs (~7uM total strand concentration) in lOmM sodium phosphate 
buffer (pH 7) at 260 nm for ODN-1 (open square), ODN-2 (filled 
square), ODN-3 (open triangle) and ODN-4 (filled triangle), and at 
275 nm for ODN-5 (open circle) and ODN-6 (filled circle). 




30 40 50 60 

Temperature (°C) 

Figure 8. PPC curves in 10 mM sodium phosphate buffer at pH 7.0: 
ODN-3 (open triangle, at 660 uM) and ODN-4 (filled triangle, at 
500 uM); dotted line and dashed lines are slopes for post-transition 
baseline for ODN-3 and ODN-4, respectively. 



Table 2. PPC results for the unfolding of ODN-3 and ODN-4" 



Pre-transition Post-transition 



ODN a(xl0 3 ) 3a/6T(xl0 5 ) a(xl0 3 ) 3a/6T(xl0 5 ) AK.dv (x 10 3 ) # K (xl0~ 3 ) AFu(cm 3 /mol) AC P (cal/K.mol) 



ODN-3 1.6 -8.1 1.0 -2.8 -8.0 3.5 -28 550 

ODN-4 1.8 -10.6 1.0 -2.7 -3.7 3.5 -13 20 



"These parameters are measured in 10 mM sodium phosphate buffer (pH 7.0) using a 580 uM solution of the ODN. The experimental errors are: AV 
(±3cm 3 /mol) and AC P (±10cal). 
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-28cm 3 /mol (ODN-3) and -13 cm 3 /mol (ODN-4), 
revealing volume contractions of the system as the tem- 
perature is increased. If the existence of voids is considered 
negligible, these volume contractions reflect changes in the 
hydration state of the ODNs (50,51), i.e. ODN unfolding 
is accompanied by a decrease in its volume of hydration. 

Furthermore, we have performed DSC melts at several 
salt concentrations to measure indirectly the heat capacity 
(AC P ) in the unfolding of each duplex. In the plots of 
Ai/cai versus T M for each ODN shown in Figure 9, the 
slopes of the resulting lines correspond to AC p of 
550cal/K.mol (ODN-3) and 20cal/K.mol (ODN-4). 
Positive heat capacity values imply exposure of hydropho- 
bic groups to the solvent upon duplex unfolding, which is 
expected due to base unstacking. The heat capacity effect 
is 28 times smaller for ODN-4 and is consistent with its 
hydrophilicity relative to ODN-3, as indicated earlier. The 
correlation of the volume contractions with the AC P 
measurements, and assuming similar random coil states, 
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Figure 9. The plots of the enthalpy versus T M for ODN-3 (unmodified, 
open triangle) and ODN-4 (modified, filled triangle) showing positive 
heat capacity effect for ODN-3. 



suggests that ODN-4 is less hydrated with fewer 
condensed ions, i.e. its folded state has a higher heat 
capacity value excluding more water molecules. 
However, as both duplexes have similar $j/S (Table 2), 
an alternative explanation is that ODN-4 has a higher 
percentage of electrostricted water molecules, i.e. its 
folded state interacts better with the solvent, changing 
its number of hydrating water molecules. 

NMR studies 

The NMR spectra of the unmodified (ODN-3) and 
modified (ODN-4) oligonucleotides of Figure 10 show 
the NOE connectivity of the purine Nl and pyrimidine 
N3 imino protons. The base imino protons were 
assigned based on their sequential connectivities in 
NOESY spectra, and these assignments were supported 
by their NOE cross-peaks to Watson-Crick base-paired 
amino protons (57). The sequential connectivities were 
obtained from base pairs G 3 -C 10 A 4 -T 9 -s> G 5 / 
oxoG 5 -C 8 C 6 -G 7 . For both duplexes, the 

imino-proton resonances of the terminal base pairs 
G 1- C'~ are lost by fast exchange with water. Because 
the T 11 imino resonance in modified duplex was missing, 
no NOEs arising from the T 11 imino proton were 
observed. The imino resonance from oxoG was shifted 
upfield (by 0.2 ppm), which reflected the effect of base 
pairing opposite C. 

Series of ID spectra for the exchangeable protons were 
recorded at 5, 15, 25, 35, 45, 55 and 65°C and are shown in 
Figure 11 A. The temperature dependences of the line 
widths for base pairs of unmodified and modified 
duplexes are compared in Figure 1 IB. The data show 
that around 45°C the N 2 -imino proton of the oxoG 5 *C 8 
base pair in modified ODN-4 duplex begins to broaden, 
compared to the G 5, C 8 base pair in unmodified duplex 
(ODN-3). At ~55°C, the oxoG 5 peak was broadened to 
the point of not being observable. Moreover, the T 9 *A 4 
imino proton of the flanking base pair was significantly 
broader by 25°C in the modified duplex, and baseline at 
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Figure 10. 'H NOESY spectrum showing resonances for the thymine and guanine imino protons and sequential NOE connectivity for the imino 
protons of the base pairs G 3, C 10 -* C 6, G 7 for unmodified duplex ODN-3 (left) and modified duplex ODN-4 (right). 
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Figure 11. (A) 'H-NMR of imino proton resonances as a function of temperature for unmodified duplex (left) and modified duplex (right). (B) 
Temperature dependence of line widths of the imino proton resonances of unmodified duplex (left) and modified duplex (right). 



35°C. Similarly, the G 3- C 10 imino proton broadened at an 
even lower temperature and near 30°C disappeared into 
the baseline. Imino proton of the T 11- A 2 base pair was 
undetectable compared to unmodified duplex. These ob- 
servations suggest that the oxoG 5 and T imino protons 
were in enhanced exchange with the solvent and indicated 
a localized destabilization of the oxoG 5- C 8 and T 9, A 4 
base pairs. This is consistent with reduced base pair 
stacking indicated by the reduced hyperchromicity in the 
UV melts and the lower enthalpic stabilization. 



DISCUSSION 

Based upon X-ray crystal and NMR structural studies, it 
has been suggested that DNA with an oxoG»C base pair is 
virtually identical to DNA with a Watson-Crick G»C 
base pair (30,31). Because of the perceived structural simi- 
larity, it was proposed that 8-oxoguanine-DNA 
glycosylase would need to find its substrate by probing 
the double helix for a characteristic electronic signature 
of the non-canonical oxoG deoxynucleoside (34-37). 
Theoretical calculations indicate that the cost for extru- 
sion of oxoG from the helix would be ~7 kcal/mol lower 
than for an unmodified G in the presence of the 



glycosylase (34). Our results demonstrate changes in free 
energy with a A AG ranging from 2 to 8 kcal/mol for un- 
modified versus oxoG modified DNA in the absence of 
protein. Thus, the local destabilization of DNA by the 
oxoG adduct may afford a thermodynamic signature 
that is exploited by repair glycosylases in the initial 
screening of the genome, a suggestion originally 
proposed by Plum and Breslauer (59). This mode of 
initial recognition suffers from being non-selective in 
that any base modification that reduces local stability 
could be 'recognized' by a DNA glycosylase. However, 
discrimination is achieved by steric and electronic con- 
straints that allows the modified base to fit into and 
be stabilized in the active site with a conformation that 
facilitates its enzymatic cleavage from the DNA 
backbone. This hypothesis is consistent with the 
proposed mode of action of UNG, which removes uracil 
from DNA (38-41). 

The change in the free energy of DNA with an oxoG»C 
base pair is accompanied by a major reduction in 
enthalpic stabilization. To obtain further insight into the 
reasons for the large change in the enthalpy term, studies 
were initiated to understand how the oxoG lesion affected 
hydration and cation association. The profound reduction 
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Table 3. Differential thermodynamic profiles for G -> oxoG and c 7 G — >• G substitutions in a DNA duplex" 



ODN 


AAG° 10 


AA// Cal 


A{TAS) 


AAw Na+ 


AAn w 


AAV 




(kcal/mol) 


(kcal/mol) 


(kcal/mol) 


(mol Na + /mol) 


(mol w/mol) 


(cm 3 /mol) 


2 versus 1 


+3.3 


+35.7 


+32.4 


+0.5 


+ 15 


n.d. 


4 versus 3 


+3.7 


+38.9 


+35.2 


+2.4 


+27 


-15 



"Values derived by subtracting the thermodynamic parameters of ODN-2 versus ODN-1 and ODN-3 versus ODN-4 in lOmM sodium phosphate 
buffer (Table 1). 



in the level of hydration associated with the oxoG»C base 
pair correlates with the reduction in enthalpic stabiliza- 
tion, and raises the question of what types of water are 
affected by oxoG. There are waters associated with the 
ionized phosphate backbone, often referred to as 
electrostricted water. These waters exchange rapidly and 
are seldom seen in high-resolution structures of DNA 
(60,61). There is no reason to suspect that oxidation of 
guanine will translate to any change in the hydration of 
the phosphate backbone since the unmodified and oxoG 
modified DNAs all show a B-DNA conformation in the 
CD. Therefore, water associated with the sugar-phos- 
phate backbone is probably not involved in the change 
in hydration and enthalpic destabilization induced by 
oxoG. Another type of water correlated with duplex for- 
mation is termed structural water that interacts with the 
polar groups in the DNA grooves. This type of water, 
which exchanges slower than electrostricted water and is 
observable in many crystal structures (62-64) forms a 
primary layer that is H-bonded with atoms that line the 
floor of the major and minor grooves. A third type of 
water associates with the hydrophobic surfaces of DNA 
that are exposed to solvent (65). Because the structure 
around each type of water varies, they have different 
partial molal volumes. 

Each base pair in unmodified B-DNA has ~20 waters 
that are part of a hydration network that stabilizes the 
duplex (66-68). This is considered structural water. 
Upon denaturation, a net release of four waters/base 
pair occurs despite the opening of the helix that exposes 
numerous polar H-binding sites for water molecules, 
which were tied up in base pairing in the folded state. 
Clearly, the bridging array of waters that line the floor 
of the grooves provides enthalpic stabilization of the 
duplex. When a modified base that affects normal 
pairing is present, the equilibrium between stacked and 
unstacked bases will be shifted to the latter, and stabilizing 
structural water will be replaced with hydrophobic water 
on the surface of the exposed planar nucleobase. Thus, the 
exposure of hydrophobic bases has an enthalpic penalty 
due to reduction in H-bonding that occurs in the first-shell 
of water molecules that line a planar aromatic surface. 
Accordingly, the number of hydrophobic water molecules 
associated with DNA will be minimized. This will result in 
an overall reduced hydration and would partially account 
for the unfavorable enthalpy for the oxoG»C base pair 
and the concomitant increase in entropy due to the net 
release of groove associated water molecules. 

We observed that the volume contractions (negative 
AV V ), measured in the PPC experiments, contrast with 



the water release measured using the osmotic stress tech- 
nique (positive A« w ). The reason for this discrepancy 
results from the complex hydration effects accompanying 
the unfolding of a nucleic acid duplex summarized above. 
In the duplex-coil transition, the net effect is the interplay 
of electrostricted, structural and hydrophobic water that 
depends on the sequence and temperature. Furthermore, 
the PPC and osmotic stress methods measure different 
types of water; the osmotic stress technique measures the 
thermodynamic number of immobilized/released struc- 
tural waters while PPC measures the effects of hydropho- 
bic, structural and electrostricted water. 

To determine the change in the type of water associated 
with DNA that occurs due to the G (ODN-3) ->■ oxoG 
(ODN-4) substitution, we created a Hess cycle with the 
two duplex-unfolding reactions and assumed that the 
random coil states at high temperatures are similar. This 
model yields the following reaction: unmodified duplex — > 
modified duplex, and the associated differential thermo- 
dynamic profiles in lOmM salt are shown in Table 3. The 
positive signs of the A AG (+3.7 kcal/mol) term and of 
the AAH— A(TAS) compensation (+35.2 kcal/mol), are 
opposite to the sign of the differential unfolding volume 
term (AAFtj = — 15cm 3 /mol). This indicates the partici- 
pation of structural water molecules, most likely a conver- 
sion of structural to electrostricted water, as determined 
previously for other processes (69-71). Furthermore, we 
also measured removal of structural water molecules 
(AA77 W = 27) and counterions (AA« Na+ = 2.4) resulting 
from the G — >■ oxoG substitution. The negative sign of 
the AAFij term indicates an increase immobilization of 
water molecules. Since both ODN-3 and ODN-4 duplexes 
have similar molar volumes (3v = 3.5 x 10 3 cm 3 /mol), 
this net hydration change involves an exchange of struc- 
tural water for electrostricted water molecules (71). The 
measured removal of 27 structural waters contributes with 
a volume expansion of M22cm 3 /mol (i.e. 27 x 4.5, where 
4.5 is the change in the molar volume in cm 3 /mol of struc- 
tural water) (72). We speculate that an immobilization of 
~43 electrostricted water molecules [(— 15+ 122cm 3 )/2.5, 
where 2.5 in cm 3 /mol is the change in the molar volume of 
electrostricted water (53,71)]. This net hydration exchange 
involving interconversion of structural to electrostricted 
water molecules due to the inclusion of oxoG may also 
contribute to the observed differential enthalpy term. 

The temperature-dependent NMR study shows 
enhanced exchange of the oxoG base with solvent and 
confirms the thermodynamic instability induced by the 
lesion and that its effect extends several base pairs in the 
3'-direction. It is not clear why the previous 
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temperature-dependent NMR structure of oxoG in the 
5'-CGC-oxoG-AATTCGCG sequence did not show a 
sequence-dependent enhanced exchange in the imino 
proton spectrum at oxoG»C (31). It is likely that 
ODN-4 shows a selective destabilization from the oxoG 
through the 5'-terminus because of the stabilizing influ- 
ence of the G/C central core that remains relatively 
intact until higher temperatures (>45°C). In the 
Dickerson dodecamer, which was the focus of the 
previous NMR study (31), there is an A 2 T 2 core that ex- 
changes at lower temperatures (>35°C) in the unmodified 
duplex, so the entire oxoG substituted duplex appears to 
be destabilized. 

What is the origin of the destabilization induced by 
oxoG? One possibility is the hydrophilicity of the oxoG 
base mentioned, which could reduce the penalty for it to 
unstack and be exposed to solvent vs. unmodified G. 
Another possibility is that the electronic nature of oxoG 
could reduce H-bonding. However, a prior study on the 
base pairing properties of oxoG with C, A, G and T using 
the 5',3'-bis silylated deoxynucleoside derivatives in 
CDCI3 solvent, which eliminates contributions due to 
base stacking and H-bonding involving the sugar ring 
hydroxyl groups, showed that oxoG»C and G»C have 
similar binding constants (33). Thus, the intrinsic 
H-bonding properties of oxoG are not perturbed due to 
electronic changes in the base. The potential preference for 
the syn- over the 'normal' a«?;-conformation around the 
glycosidic bond was also ruled out by the NMR studies in 
CDCI3, although the barrier for the two conformations 
may be different in double-stranded DNA. However, 
neither the NMR or crystal structure studies indicate the 
presence of a .mz-conformation for oxoG when it is paired 
with dC (29-33). 

An important consequence of the change from G«C to 
oxoG»C is the loss of a major groove cation-binding site 
due to the replacement of an electronegative N7 atom with 
an electropositive amide-like N-H. Monovalent cations in 
high-resolution crystal structures are generally observed 
near the major groove edge of G«C base pairs close to 
either the N7 and 06 atoms or equi-distant between 
the two (73,74). Another high-occupancy monovalent 
cation binding site is midway between stacked interstrand 
G's at 5'-GC sequences (69,70), which is the environ- 
ment of oxoG in ODN-4. We have previously reported 
that 7-deazaguanine (c 7 G) (i.e. Figure 1, ODN-7 
and ODN-8) causes a destabilization of 1.3-2.5 kcal/mol 
(over a range of 10-200 mM NaCl) that involves reduced 
enthalpic stabilization and accompanying reduction in hy- 
dration and cation binding (42,75). The destabilization did 
not occur at the c 7 G»C base pair but at the neighboring 
normal C»G base pair. It was shown that the cation, 
which buttresses the interstrand G's at the 5'-GC 
sequence, was eliminated in the crystal structure (43). 
Similar to what has been reported for oxoG»C base 
pairing (33), the c 7 G»C base pair did not impact 
H-bonding affinity. We propose that the loss of the elec- 
trostatic association of the cation with the major groove 
edge near 5'GC destabilizes the DNA, which is revealed in 
reduced base stacking and reduced hydration. 



SUMMARY 

Although DNA with oxoG lesions appears to have a 
similar structure to natural DNA at low temperature, 
the thermodynamic and temperature-dependent NMR 
studies indicate that the environment around the lesion 
is significantly altered. Changes in cation binding and hy- 
dration accompany the enthalpic destabilization of DNA 
by the presence of oxoG. 
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